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EFFECTS OF WING ELASTICITY ON THE AERODYNAMIC 
CHARACTERISTICS OF AN AIRPLANE CONFIGURATION HAVING 
45° SWEPTBACK WINGS AS OBTAINED FROM FREE-FLIGHT 
ROCKET-MODEL TESTS AT TRANSONIC SPEEDS 
By A. James Vitale 


SUMMARY 


Flight tests at transonic speeds have been made of tvo rocket- 
propelled models of an airplane configuration having k5^ svepthack wings 
of aspect ratio h and with different stiffnesses. The longitudinal sta- 
bility, lift, and drag characteristics were determined by analysis of 
the response of the models to disturbances in pitch. A comparison was 
made between experimental and calculated effects of aeroelasticity on 
lift and longitudinal stability. 

A straight-line extrapolation of the lift-curve slopes of the two 
wings of different stiffness to obtain the rigid-wing lift-curve slope 
was found to be satisfactory in the range of wing stiffness tested. The 
loss in lift-curve slope of the more flexible wing found experimentally 
agreed fairly well with that calculated by a modified strip-theory method. 
A combination of the effects of a high-tail position and an unstable 
pitching-moment curve of the wing at high lift coefficients and tran- 
sonic speeds caused both models to execute a violent pitch-up maneuver. 

The induced drag parameter dCp/dCL^ was found to be nonlinear with 
lift coefficient and the model having the more flexible wing had lower 
drag due to lift over the Mach number range of 0.97 to I. 30 . 


INTRODICTION 


The results presented in this paper are part of a general research 
program utilizing rocket -propelled models to investigate the effect of 
wing plan form on transonic and supersonic longitudinal stability, lift, 
and drag characteristics. Previous models tested have included several 
thin straight wings and triangular wings (refs. 1 to 6 ) . All the models 
tested have had similar fuselage and tail surfaces. The wing-off charac- 
teristics of the models tested are presented in reference 2 . 


V 


CONFIDENTIAL 


2 


CONFIDENTIAL 


NACA RM L52L30 


In this investigation tvo models were tested which had wings of 
different stiffnesses. The two wings differed in stiffness hy reason 
of the material of construction, one having been constructed of solid 
steel nnrl the other of solid duralumin. Since the only difference 
between the two models was wing stiffness, differences in the test 
results may be attributed to aeroelastic effects. The wings which had 
an aspect ratio of U and a taper ratio of 0.60 were swept back k5° at 
the quarter-chord line. The wing airfoil sections parallel to free 
stream were NACA 65A006. 

An all-movable horizontal tail was used for longitudinal control 
on the models. During the flight the horizontal tall was deflected in 
an approximate sqijare-wave program between the limits of 0.l4° and - 3 . 0 ° 
for the steel-wing model and - 0 . 72 ° and -k. 60 ° for the duralumin-wing 
model. The pulsed-control technique discussed in detailed form in refer- 
ence 1 was used to obtain and analyze the data. 

In addition to experimental results, the losses in lift and longi- 
tudinal stability due to aeroelastic deflection were calculated for the 
duralumin-wing model by the method outlined in the appendix of this 
report . 

The models were flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 


Cn 

normal-force coefficient. 

an W/S 
g q 

Cc 

chord-force coefficient. 

Ifl w/s 
g q. 

<^L 

lift coefficient. 

Cn cos 

cx ”• Cq 1 

CD 

drag coefficient. 

Cc COB 

a + Cjj 1 

Cjn 

pitching-moment coefficient 

Kw 

wing normal force. 

lb 



Nv 


^Ny *^Lw wing lift coefficient, 


^n 


qS 

normal acceleration as obtained from accelerometer, 
ft/sec/sec 


CONFIDENTIAL 


NACA RM L52L30 


CONFrOENTIAL 


3 




W 

g 

P 

Po 


<1 

M 

R 

S 

5 

c 

b/2 

y 

y' 

L 

e 

(9/L)^ef 


a 

5 

7 

A 

P 


longitudinal acceleration as obtained from accelerometer, 
ft/sec/sec 

model veight, lb 

acceleration of gravity, ft/sec/sec 

free-stream static pressure, Ib/sq ft 

standard sea-level static pressure, 2116 lb/sq_ ft 

dynamic pressure, ^-pM^ 

Mach number 

Reynolds nimiber, based on wing mean aerodynamic chord 

wing area, including the area enclosed within the 
fuselage, sq ft 

wing mean aerodynamic chord, ft 
wing chord, ft 
wing semlspan, ft 

lateral distance from fuselage center line, ft 
lateral distance from fuselage side, ft 
load applied, lb 

local wing twist angle produced by L, radians; also 
angle of pitch, rad.ians 

etructxiral influence coefficient at spanwise center of 
pressure 

angle of attack, deg 
horizontal-tail deflection, deg 
flight-path angle, deg 
aspect ratio 

period of oscillations, sec 
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'^1/2 


time to damp to one-half amplitude, sec 


I 


distance from center of gravity to nose normal 
accelerometer, ft 


Subscripts : 
d = — - ^ 

57.3 

= ^ 3 - 

dt 2V 


w 

wing 

e 

elastic 

r 

rigid 

cp 

center of pressure 


MODELS AND APPARATUS 
Models 


A three-view drawing of the models used in this investigation is 
shown in figure 1. Photographs of the models are given in figure 2. 

The fuselage has a cylindrical center section with identical ogival 
type of nose and tail sections defined hy the ordinates given in table I 
The body fineness ratio is I 3 .O. 

The two models tested had identical wings except for material of 
construction. One model had a wing constructed of solid steel, and the 
other model had a wing constructed of solid duralumin. The wing airfoil 
sections parallel to free stream were NACA 6SA006 airfoil sections. 

The horizontal tail was constructed of solid duralumin for both 
models and was mounted on a ball bearing built into the vertical tail 
at a height of 0.508 semispan above the wing-chord plane extended. A 
hydraulic power system was used to deflect the horizontal tail in an 
approximate square-wave program. The gap between the vertical tail and 
the root of the horizontal tail was sealed by means of a wiper-type seal 
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The mass characteristics of the two models tested are given in 
the following table: 



Steel- wing model 

Dural- wing model 

w 

ly 

c.g. position 

lk6 pounds 
10.30 slug-ft 2 
2k .8 percent MAC 

119 pounds 
10.95 slug-ft2 

k.6 percent MAC 


Instrimientation 

The models were equipped with telemeters transmitting measurements 
of normal acceleration at the model center of gravity, angle of attack, 
longitudinal acceleration, wing normal force, control position, total 
pressinre, reference static pressure, and transverse acceleration. The 
duralumin-wing model had a normal accelerometer located 2.68 feet ahead 
of the model center of gravity in addition to the center-of-gravity 
normal accelerometer. This arrangement permitted the determination of 
instantaneous pitching moment for the duralimiin-wing model. 

A vane-type instrument mounted on a sting extending frqm the nose 
of the model (fig. 2) was used to measure angle of attack. For the 
'duralumin-wing model this sting was deflected downward 5 ° to permit 
positive angle -of -attack measurements up to 20°. The angle-of -attack 
indicator is described in more detail in reference 7 * 

The wing was mounted on a beam-type balance calibrated to give 
measurements of wing normal force. The balance is described in more 
detail in reference 3 * addition to the inductance -gage pickup 
described in reference 3 ^ the balance for the duralumin-wing model 
was equipped with a strain-gage pickup for development purposes. 

The total-pressure pickup was mounted on a small strut below the 
fuselage. The static -pres sure orifice was located k.9 inches behind 
the beginning of the cylindrical portion of the fuselage. A calibra- 
tion of the reference static pressure for an angle of attack of 0 ° was 
obtained from previous instrumentation models. 

Atmospheric conditions were determined from radiosondes released 
shortly after the flights . Fixed and manually operated l 6 -millimeter 
motion-picture cameras were used to photograph the laiinchings and first 
portions of the flights. Additional ground equipment consisted of a 
CW Doppler radar unit for obtaining model velocity, a tracking radar 
unit for obtaining model position in space, and special (spinsonde) 
radio equipment for obtaining model rolling velocity. 
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TESTS 

Preflight Measxirements 


Prior to flight testing the structural influence coefficients were 
obtained for the duralumin wing. These coefficients are shown in fig- 
ure 3 as the twist in the free-stream direction per \init load applied 
along the 25- and 50-percent chord lines • 

With instruments installed, the models were suspended hy shock 
cords and vibrated in the pitch plane by an electromagnetic shaker and 
also by striking the wing and fuselage. The following model natural 
frequencies and modes of vibration were determined from the telemeter 
records taken during these ground tests and from visual observations of 
the model vibrating. 


Wing mode 

Frequency, cps 

Steel wing 

Duralumin wing 

First bending 

kk 

46 

Second bending 

212 

206 

Torsion 

340 

350 


Model .Tests 

The models were lamiched at approximately 60° from the horizontal 
by means of a mobile launcher shown in figure L. A 6- inch-diameter 
solid-fuel ABL Deacon rocket motor was used to boost each model to 
maximum velocity. The models contained no sustaining rocket motors 
and experienced decelerating flight after separation from their boosters. 
During the decelerating portion of the flights, the models experienced 
short-period oscillations in angle of attack, normal acceleration, longi- 
tudinal acceleration, and wing normal force following each horizontal- 
tail deflection. 

The measured angle of attack was corrected to the angle of attack 
at the model center of gravity by the method of reference 7- The wing 
normal-force measTurements were corrected for inertia effects to give 
aerodynamic forces. 

The Mach nvunbers and dynamic pressures were calculated from telem- 
etered total pressure and free-stream static pressure obtained from the 
combination of radiosonde and tracking radar data. The velocity obtained 
from the CW Doppler radar unit furnished a check at peak velocity on the 
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Mach numbers obtained from the pressure data. For use in comparing the 
aeroelastic properties of the wings in the present investigation with 
other results the values of free-stream pressure obtained during the 
model flights divided by standard sea-level pressure are presented in 
figure 5 as a function of Mach number. 

Reynolds numbers (based on wing mean aerodynamic chord for each 
model) obtained during the flights are shown in figure 6. 


Time histories of the model short-period oscillations in angle of 
attack, normal acceleration, longitudinal acceleration, and wing normal 
force were analyzed by the method of reference 1 to obtain the longi- 
tudinal stability, lift, and drag characteristics for each model. 

In addition to obtaining longitudinal stability characteristics by 
the method of reference 1 the instantaneous pitching moment was measured 
for the d\iralumin-wing model by means of two normal accelerometers located 
at the center of gravity and in the nose of the model. The difference in 
the two normal-accelerometer readings is proportional to the angular 
acceleration in pitch by the relation: 


The pitching moment due to angle of attack is then calculated by the 
following equation: 


For the steel-wing model, a nose accelerometer was not used but the 
instantaneous pitching moments were calculated from the flight data by 
use of .a differentiation process. The angular acceleration in pitch was 
obtained from the following equation: 


The quantity d was obtained by differentiating the measured a curve 
and the quantity y was calculated from the measured accelerations at 
the model center of gravity, the gravity component being neglected. 


Analysis 


^%iose ~ ^’’cg 


I 



dt dt 
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This procedure is the same procedure used in applying the corrections 
to the angle of attack measured at some point other than the center of 
gravity (ref. 7 ). Neglecting the gravity component has a negligible 
effect for these models. 

The effects of aeroelastic distortion of the duralumin wing on 
lift and stability characteristics were calculated by means of a modi- 
fied strip theory and the structural influence coefficients shown in 
figure 3. The method of calculating the aeroelastic effects is out- 
lined in the appendix of this report. 

Data were obtained for the steel-wing model from Mach numbers 0.97 
to 1.27 and for the duralumin -wing model from Mach numbers O.92 to 1 . 33 . 
At a Mach number of 0.97 for the steel-wing model and 0.92 for the 
duralumin-wing model, the models experienced a pitch-up maneuver to 
high angles of attack and were not able to recover; thus further data 
could not be obtained. 


Accuracy 

The following tables give estimated values of the possible systematic 
errors in the absolute values of Cl, Cp, and ns affected by the 

accelerometer and wing-balance calibration ranges ; 

For the steel-wing model. 


M 

^L 



0.97 

1.15 

1.30 

±0.020 
±.013 
+ .009 

±0.016 

±.010 

±.008 

±o.ook 

±.003 

±.002 


and for the dural\jmin-wing model, 


M 

Z^L 


^D 

0.92 

1.15 

1.35 

± 0.015 

±.009 

±.006 

±0.011 

±.006 

+ .005 

tO.OOk 

±.0023 

±.0016 
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Further errors in aerodynamic coefficients may arise because of 
possible dynamic -pressure inaccuracies which are approximately twice 
as great as the errors in Mach number. The Mach numbers are estimated 
to be accurate to ±1 percent at supersonic speeds and ±2 percent at 
subsonic speeds. 

The exposed-wing lift data indicate a further error in the absolute 
value of the wing normal-force measurements for the duralumin-wing model. 
As stated in the instrumentation section, the wing balance for the 
duralumin -wing model was equipped with a strain-gage pickup for develop- 
ment purposes . The data from the strain-gage pickup do not indicate 
the zero-shift of about 200 pounds shown in figure 8 Tor the exposed 
duralumin-wing lift data from the inductance -gage pickup. This zero 
shift does not affect the incremental values or slopes which eire, in 
general, more accurate than the absolute values. 

The errors in the measured angles of attack and elevator deflec- 
tions should not vary with Mach number because they are not dependent 
on dynamic pressure. The greatest possible error in angle of attack is 
caused by possible aerodynamic asymmetry of the angle-of -attack vane 
which is not detectable prior to flight. The following table gives 
estimated values for the errors in absolute values of a and 6 as 
affected by instrument calibration ranges: 



Ax 

A6 

Steel-^ing model 
Duralumin-wing model 

+ 0 . 17 ° 

+ .25° 

+0.1° 
+ .1° 


RESULTS AND DISCUSSION 
Lift 


The basic lift data are shown in figures 7 and 8 for the duralumin- 
and steel-wing models at several Mach mmibers and horizontal tail deflec- 
tions. In figure 7> the total lift coefficients are plotted against 
angle of attack for both models . It is seen from these curves that the 
lift-curve slope is nonlinear with lift coefficient. In figure 7(a) it 
is seen that the lift-curve slope near zero lift at Mach numbers 1.02 
and 1.22 is about 15 percent lower than the lift-curve slope near lift 
coefficient of O.3O. At higher angles of attack as shown by the lift 
curves for Mach numbers 0.97, 1.12, and 1.27, the lift-curve slopes 
again decrease due to separation over the wing tips. Shown in figure 7("b) 
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is the lift curve for the duralumin-wing model at a Mach number of O .92 
and over sui angle -of -attack range from 4° to l8°. The lift coeffi- 
cient 0.70 at which the pitch-up occvirred is well below the maximum 
lift coefficient obtainable at a Mach number of O. 92 . 

In figure 8 the exposed-wing lift coefficients are plotted against 
angle of attack for both models. The wing balance measures normal force 
only; however, the difference between normal force and lift force was 
estimated to be small for the exposed wing and the data are plotted as 
lift coefficients. The same nonlinear trends described in the previous 
section on total airplane lift are seen to apply also to the exposed- 
wing lift. This phenomenon of increase in slope at 0.20 to 0.40 lift 
coefficient and then decrease in slope at higher lift coefficients was 
also found in wind-tunnel tests on this wing (refs. 8 and 9 ). 

The variation of lift-curve slope with Mach number is shown in 
figure 9 for the exposed wings and complete models. The slopes are 
taken near a lift coefficient of O. 3 O since data were available at this 
lift coefficient for both horizontal-tail deflections. The duralumin- 
wing lift-curve slopes for the complete model and exposed wing increase 
rapidly from a Mach number of 0.92 to 0.95 and then gradually decrease 
as the Mach number increases. The steel-wing lift-curve slopes show the 
same variation with Mach number but, as would be expected, have higher 
values than for the more flexible duralumin wing. Also shown in fig- 
ure 9 is the lift-curve slope of the exposed rigid wing as obtained by 
the method described in the section entitled "Aeroelastic Effects on 
Lift . " 


The ratio of the exposed-wing lift to complete-model lift is shown 
in figure 10. The values were obtained by dividing the exposed -wing 
lift-curve slopes by the complete-model lift-curve slopes shown in fig- 
ure 9 for a lift coefficient of O. 3 O. The contribution of the wing to 
the complete -model lift is approximately constant over the Mach number 
range covered for both models. The ratio of exposed-wing area to total 
wing area is 0.78 for the models tested. For the steel-wing model the 
ratio of exposed-wing lift to total lift varies from 0.75 at a Mach num- 
ber of 1.0 to 0.71 at a Mach number of 1.27 . When the tail contribution 
to the total lift is taken into account the results for the steel wing 
show that the usual assumption that the fuselage carries that part of 
the load represented by the wing area intercepted by the fuselage gives 
nearly correct answers over the Mach number range 0.97 to 1.27. The 
exposed duralumin wing carried a slightly smaller percentage of the 
total lift than the exposed steel wing from Mach numbers 1.0 to I. 3 O. 
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Model Vitrations 

For the duralumin-wing model there was a high-frequency oscillation 
of 180 cycles per second superimposed on the telemetered traces of angle 
of attack and nose normal acceleration from Mach numbers 1A6 to 1.25. 

This type of nose vibration which started near the time of booster -rocket 
burnout also occurred on a similar model having a considerably different 
wing tested in reference 6. 

Near the beginning of the pitch-up maneuver of the steel-wing model 
(discussed in the static stability section) at a Mach number of 0.9^5 
and lift coefficient of 0.64 the wing normal-force and normal-acceleration 
measurements indicated a vibration corresponding to the wing first-bending 
frequency which was probably caused by wing buffeting. A similar vibra- 
tion occurred at the beginning of the pitch-up maneuver of the duralumin- 
wing model at a Mach number of 0-92 and lift coefficient of O.7O. At a 
Mach number of 0.955 and lift coefficient of 0.56, there also was an 
indication of wing buffeting for the duralumin-wing model. 


Aeroelastic Effects on Lift 

For use in calculating the loss in lift of the duralumin wing, the 
structural influence coefficients were measured before the model was 
tested. These influence coefficients axe shown in figiure 3 concen- 
trated loads along the 0.25 chord line and the 0.50 chord line. It was 
assumed that structural influence coefficients of the steel wing were 
related to the values measured for the duralumin wing by the ratio of 
the modulus of elasticity of duralumin to that of steel. One method of 
using these structural influence coefficients is shown in figure 11. 

Data obtained from reference 10 at a Mach number of 1.11 for a rocket 
model equipped with the same wing constructed of wood with duralumin 
inlays is compared with the complete-model lift-curve slopes of the 
steel and duralumin -wing models. The abscissa of figure 11 is the 
parameter qS(0/L)pg^ where is the structural influence 

coefficient from figure 3 the spanwise center of pressure of the 
exposed rigid wing for the load also applied at this spanwise location. 
Since the exposed-wing lift-curve slope was not measured in the test 
of reference 10 for the wood-dufalumin wing and the ratio of tail area 
to wing axea is nearly the same as for the steel- and duralumin-wing 
models, the complete-model lift-curve slopes are plotted against the 
parameter qS(a/L)pgf . When influence coefficients for both the 

0.25-chord loading and the 0.50-chord loading are used, figure 11 shows 
that a straight-line extrapolation is possible over the range of stiff- 
nesses for the three wings tested to obtain the value of the rigid lift- 
curve slope. The same method was used in figure 12 to obtain the rigid 
exposed -wing lift-curve slope at several Mach numbers. The variation of 
the rigid exposed-wing lift-curve slope with Mach number is shown in 
figure 9. 
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The loss in lift of the dviralumin wing was calculated hy means of 
the modified-str ip-theory method outlined in the appendix. The method 
used allows the exposed wing to be divided into any number of sections 
and also allows any desired type of span loading to be used. Shown in 
figure 13 is the ratio of the elastic- to rigid-exposed -wing lift-curve 
slope of the duralumin wing calculated for two -point and five -point span 
loadings from a span— loading curve of reference 11. Since the exact 
location of the chordwise center of pressure was not known, the calcu- 
lations were made by using the influence coefficients for both the 
0.25-chord loading and the 0.50-chord loading. Pressure measurements 
on a similar wing tested in reference 12 show that these two positions 
should bracket the chordwise center-of -pressure variation along the 
span and that for this wing the 0.50-chord loading is more correct over 
the Mach number range 1.0 to I.30. The results of figure I3 show that, 
for the wing tested, there is little difference in the ratio of the 
elastic to rigid lift-curve slope calculated by a two-point loading and 
by a five-point loading. The results from a calculation using a two- 
point loading from a trapezoidal span-load curve indicated only a 
1-percent greater lift loss than a two-point solution using a span- 
load curve of reference 11. 

The effect of wing-inertia loading which acts in the opposite 
direction to aerodynamic loading was included in the lift-loss calcula- 
tions shown in figure l4. These calculations were made using equation 7 
of the appendix for a two-point loading which includes the wing-inertia 
effect. The weight of one panel of the exposed duralumin wing was 
5.93 pounds and the effect of including wing inertia in a two-point 
solution reduced the calculated lift loss by about 2 percent in the 
Mach number range 0.97 to 1.27 where the elastic -exposed -wing lift- 
curve slope varies from 90 to 80 percent of the rigid value. 

A comparison between the calculated values and experimental values 
of the ratio of the exposed-dural\jmin-wing and exposed-rigid-wing lift- 
curve slopes of figure 9 is also shown in figure lU. At a Mach number 
of 0.97 the experimental ratio agrees with the calculated value for a 
0.25-chord loading, and over the Mach number range 1.02 to 1.27 the 
experimental ratios agree with the calculated values for a 0.50-chord 
loading. 

The loss in lift for the duralumin wing was also calculated from 
the charts of reference I3 at subsonic speeds and the values obtained 
agree fairly well with the calculations made by the method of the 
appendix at these same speeds. The charts were not used for higher 
speeds since reference I3 states that the charts are less reliable when 
the component of Mach niimber perpendicular to the leading edge of the 
wing is transonic. 
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The wind-tunnel investigation of reference 14 for a wing-fuselage 
combination having steel and duralumin wings similar to those of this 
report showed less reduction in lift for the duralumin-wing — fuselage 
combination than shown in figure l4 for the exposed duralumin wing of 
this investigation. The dynamic pressures for the investigation of 
reference l4^ however, varied from only 400 to 850 pounds' per square 
foot as compared with 1100 to 2500 pounds per sqiiare foot for the 
rocket-model investigation. 


Static Stability 

The basic pitching-moment data are shown in figttre 15 for the steel- 
wing model at a center -of-gravity position of 24.8 percent of the mean 
aerodynamic chord, and in figure l6 for the duralumin-wing model at a 
center-of-gravity position of 4.60 percent of the mean aerodynamic chord. 
The pitching -moment coefficients are plotted against lift coefficient 
over the first one and one-half cycles of the model short-period oscil- 
lation following each horizontal-tail deflection. The measured pitching 
moments have been corrected for the damping-in-pitch derivatives, and 
the data of figures 15 and l6 do not show any consistent difference for 
increasing and decreasing angle of attack. As might be expected, the 
data obtained by using the double differentiation procedure show more 
scatter than those obtained from the double accelerometer method. The 
pitching-moment data at a Mach number of O.97 and horizontal-tail deflec- 
tion of -4. 60° in figirre 16 include the effect of a trim change at this 
Mach number. The change in trim lift coefficient is from 0.46 at a Mach 
number of 0.975 to 0.52 at a Mach number of O.965. 

At a Mach number of 0.97 the steel-wing model experienced a pitch-up 
maneuver to high angles of attack which was so violent that the model was 
not able to recover. At the peak of the maneuver the model developed a 
high rate of roll of about 15 radians per second; this high rate indi- 
cated that one wing stalled befbre the other wing. The pitching -moment 
curve (fig. 15) for this part of the flight could not be calculated above 
a lift coefficient of 0.82 because the angle -of -at tack vane hit a stop 
at about 12°, although the model continued to pitch up to considerably 
higher angles. This pitch-up corresponds to the unstable break in the 
pitching -moment curve for this wing alone at high lift coefficients and 
transonic speeds. With the tail placed in the high position of 0.508 
semispan above the wing-chord plane extended, the tail contribution to 
the longitudinal stability was not sufficient to overcome the imstable 
pitching -moment-curve §lope of the wing at high lift coefficients. 

Further discussion of the effect of high tail position is found in 
reference 15. 

The duralumin-wing model was tested with a center-of-gravity posi- 
tion approximately 20 percent of the mean aerodynamic chord more forward 
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than the steel-wing model. The pitching-moment data of figure l6 indi- 
cate that, at a Mach number of 0.97 and lift coefficient of O.73, the 
dinralumin-wing model approached the unstable break in the cui*ve but did 
not overshoot into the unstable region. At a Mach n\jmber of O.92, how- 
ever, the duralumin-wing model experienced the same pitch-up maneuver 
described previously for the steel-wing model. The pitching -moment 
curve for this maneuver (fig. 16(a)) has a very unstable break at a 
lift coefficient of O.JO and then indicates stable slope at 0.85 lift 
coefficient. The pitching-moment curve for the dural\jmin-wing model at 
a Mach number of 0.97 was converted to the center-of -gravity position 
for the steel-wing model and is shown by a dashed line in figure 15(a) 
for comparison with the data on the steel-wing model at the same Mach 
number. The curves are similar in shape and the displacement of the 
curves from each other is accounted for by the difference in horizontal- 
tail deflection. 

The measured periods of oscillation of the angle of attack are 
shown in figure 17 . The data shown for the duralumin-wing model are 
for a center-of -gravity position of 4.6 percent of the mean aerodynamic 
chord, and for the steel-wing model at 24.8 percent of the mean aero- 
dynamic chord. The periods of oscillation were converted to the longi- 
tudinal stability derivative Cmo, by the method described in reference 1. 
Since the two models were tested with different center-of -gravity posi- 
tions, the static stability is compared in figure I8 on the basis of 
aerodynamic -center position. The duralumin-wing model shows a rapid 
rearward movement in aerodynamic center from 45 percent of the mean 
aerodynamic chord at a Mach mmiber of O.93 to 60 percent at a Mach num- 
ber of 0.97. From a Mach number of 0.97 to I.33 the aerodynamic center 
gradually moved rearward from 60 percent to 72 percent of the mean aero- 
dynamic chord. Aerodynamic -center positions obtained from the pitching- 
moment data of figures 15 and 16 near trim conditions for the two models 
agree very well with the values calculated from the period of oscilla- 
tions except for one point representing a ciirve obtained from the double 
differentiation procedure. 

The aerodynamic -center position for the steel-wing model has the 
same variation with Mach number as the duralumin-wing model over the 
Mach number range 1.02 to 1.2?. Over this Mach number range, however, 
the aerodynamic -center position is 3 to 6 percent more rearward for the 
steel-wing model. From the strip-theory method of the appendix the for- 
ward movement in aerodynsunic -center position due to the inboard movement 
of the span load was calculated for the duralumin-wing model. A com- 
parison between the calculation made for a 0.50-chord loading and the 
measured difference in aerodynamic -center position for the steel- and 
duralumin-wing models in figure 18 shows that nearly all of the loss in 
stability of the duralumin-wing model may be accounted for by the inboard 
movement of load on the wing. 
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Damping in Pitch 

The time required for the pitch oscillations to damp to one-half 
amplitude is shovm in figure 19(a) for the steel-wing model and in fig- 
irre 19(h) for the duralumin-wing model. The data converted to the 
damping-in-pitch derivative Cm^ + Cm^ are shown in figure 19(c ) . The 

duralumin -wing -model data of figure 19 are presented for a center-of- 
gravity position 20 percent of the mean aerodynamic chord more forward 
than for the steel-wing model. The difference in center-bf-gravity 
position was calculated to have a small effect on the damping-in-pitch 
derivatives of figure 19(c). Also shown in figure 19(c) is the damping- 
in-pltch derivative from the wing-off model of reference 2 converted to 
the dimensions and average center -of-gravity position of the steel- and 
dviralumin-wing models. 

The pitch-damping derivative for the duralumin-wing model shows a 
rapid variation with Mach number in the transonic speed range. This 
same variation would probably have occurred for the steel-wing model 
also, but damping data were not obtained for the steel-wing model below 
a Mach number of 1.0. The lower value of C^^ + Cm^ for the duralumin- 
wing model from Mach numbers 1.02 to 1.2 may not be due to aeroelasticity 
and is more likely an indication of the accuracy of the damping- in-pitch 
derivative. Since the damping derivative for the wing-off model is 
essentially Cm^ of the' tail, the data of figiore 19(c) show that the 

wing contribution to the damping derivative Cmq + is very small 

over the Mach niomber range 1.15 to I.30. 


Drag 

The basic drag data for the steel- and duralumin -wing models is 
shown in figures 20 and 21 for several Mach numbers and horizontal-tail 
deflections. The drag data for the duralumin -wing model shown in fig- 
ure 21(c) for an average Mach number of O.97 were plotted over the 
Mach number range O.965 to O.975 and show that because of the tran- 
sonic drag rise there is a noticeable change in the level of the drag 
polars with a small change in speed. 

The minimum drag values obtained for both models are compared with 
the minimum drag of the wing-off model of reference 2 in figure 22, 

Also shown in figure 22 is the wing-plus-interference drag obtained 
from the large-scale drag test of reference I6. The wing-plus- interference 
drag obtained from the steel- and duralumin-wing-model tests compares 
favorably with the values of reference 16 at supersonic speeds. Because 
of the small size of the wing relative to the body and the high-drag body 
and tail, the wing drag represents a small percentage of the drag of the 
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steel- and duralumin-wing models. Thus, at a Mach number of 0.9^, the 
wing drag is about the same magnitude as the accuracy of the acceler- 
ometer readings of the steel- and duralumin -wing models. 

The drag due to lift of the models tested is presented in fig- 
ure 23 (a) for lift coefficients of 0.50 and in figure 23 (b) for lift 
coefficients of 0.15 since the induced drag was found to be nonlinear 
with lift coefficient. Also shown in figure 23 is the ideal induced- 
drag factor l/jtA and the drag due to lift for no leading-edge suc- 
tion 1^57.30]^. For the 1^57.3Ci^ comparison, the lift-curve slopes 

were measured near the trim-lift coefficient corresponding to elevator 
deflections of - 3 . 0 ° smd 0 .lU° for the steel-wing model and -L. 60 ° and 
- 0 . 72 ° for the duralumin-wing model. 

Near lift coefficients of 0.50 the duralumin-wing model shows a 
reduction in induced drag of about 16 percent compared to 1^57.3 CLq_ 

over the Mach number range O .92 to I. 33 . The data for the steel-wing 
model show no reduction in induced drag from Mach numbers of 1.05 to 
1.27 for lift coefficients of 0.50. 

At low lift coefficients near 0.15, the duralumin-wing model shows 
a reduction in dCp/dCL^ of about 30 percent from Mach numbers 1.02 

to 1.25 and 42 percent at a Mach number of 0.93 compared to the fac- 
tor 1 / 57 . 301 ^. The steel-wing model shows a reduction of I 8 percent 

in dCp/dCL^ at low lift coefficients over the Mach number range 1.02 

to 1.22. The wind-tunnel investigation of reference 14 for lower 
dynamic pressures than the present investigation indicates only slightly 
lower drag for the duralumin-wing — fuselage combination than for the 
steel -wing— fuselage combination. 

The lower induced drag for the duralumin-wing model suggests the 
possibility of a more efficient span loading for the flexible wing as 
a result of the inboard movement of the load when the wing deflects. 


CONCLUSIONS 


The flight tests of two rocket-propelled models having 45° swept- 
back wings of aspect ratio 4 and with different stiffnesses indicated 
the following results: • 

1. The loss in lift-curve slope due to aeroelastic distortion found 
experimentally agrees very well with that predicted by a modified-str ip- 
theory method of calculation. 
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2. In the stiffness range investigated, a straight-line extrapola- 
tion of the lift-curve slopes of wings with different stiffnesses can 
be used to obtain the value for a rigid wing. 

3. The exposed duralumin wing carried a slightly smaller percentage 
of the total airplane lift than the exposed steel wing from Mach num- 
bers 1.0 to 1.30. 

4. Because of the combined effects of an unstable break in the 
pitching-moment curve of the wing and a high tail position, the two 
airplane configurations tested experienced violent pitch-up maneuvers 
at transonic speeds and lift coefficients below the maximum. 

5. The more forward aerodynamic -center position for the airplane 
configuration having a duralumin wing was mostly due to the inboard 
movement of the spanwise center of pressure. 

6. Over the Mach number range 1.15 to I.30, the wings tested con- 
tributed very little to the damping-in-pitch derivative Cmq + Cm^. 

7* The induced -drag parameter dCj)/dCL^ was nonlinear with lift 
coefficient such that both models tested had higher values of dCp/dCL2 
at a lift coefficient of 0.50 than at a lift coefficient of 0.15. 

8. The airplane configuration having a duralumin wing had lower 
drag due to lift than the steel-wing model; this fact suggests the 
possibility of a more efficient span loading for the more flexible wing. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va. 


CONFIDENTIAL 


18 


CONFIDENTIAL 


NACA EM L52L30 


APPENDIX 

EFFECTS OF AEROELASTIC DEFLECTION 


The method used to estimate the change in lift and center of pres- 
sure due to elastic deflection of the wing is a form of strip theory. 

The information needed to apply this method is data on twist along the 
wing due to unit loads applied along the centers of pressure (taken to 
be along the 0.25 or 0.50 chord line), an assumed rigid-wing span- 
loading curve, and an estimate of the rigid-wing lift-curve slope. 

The portion of the half wing which is elastic (that portion outside 
of the fuselage) is divided into as many psirts as desired. The rigid lift 
per unit angle of attack of these sections is then estimated from the span 
loading. The increment in lift due to twist at the section is then assumed 
to be given by this lift-curve slope multiplied by the twist. Since the 
twist at a given section depends on the final equilibrium lift at all sec- 
tions, a set of n simultaneous equations (equal in number to the number 
of sections the wing is divided into) must be solved; this solution gives 
the resultant elastic lift at each section. The equations are set up in 
ratio form with the elastic lift at each section given as a ratio of the 
rigid lift at the section. The independent variable used in the equations 
is the rigid lift-curve slope times the dynamic pressure C];Q_^q. 

Sketches and equations for the exposed-wing and span-loading curve 
divided into five sections are given below. 


I 



Exposed wing panel divided into five sections. 
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^-0 ^l t ]2 Hk T]^=l 

side ^ ^ 


Span-load curve divided into five sections on exposed wing. 



side 
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where L's are the loads on each section of the exposed wing, t] = 

b/2 


K' = 


Q is defined as q, and Kq, Kg, K^, 


CZC 

*^L*^av 


dTi 


^fus . 
side 

K4, and are defined as follows: 


Ki = 


K' # 


'T)l / C7C 


side 


dr\ 


Kc: = 


K' t 


c^c 


2 /7^i^\CLCav/rigid 


dq 


> 


( 2 ) 


The ratio of the twist at each section to the rigid angle of attack 
can be derived from equations (l) and (2) for the following expressions: 


Axq Lq 

(xrKqQ 


1 




^5 _ L5 ^ 

aj* (Xj-K^Q 


(3) 


The influence coefficient, 0qj> is defined as the twist at sta- 
tion i due to a unit load applied at station j, and the equations 
for the twist of each wing section due to loads Lq, Lg, . . . L5 are 
as follows : 
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CLr^ 



^2 

ap 



L4 

if®!'* 



. ile„ t i£e,2 * %„ + iie„ i5< 

“r “r Oj. Or 23 ^2k ^ 5^ 


?T^25 




^^^5 Lp Lp LIl Lc: 

— = 55:^51 + ^52 + 5^53 + + ^55 


y 


When equations ( 3 ) are substituted into equations (4) for Axi/cxj,, 
txLzf<^} • • • Axs/ar, the following set of five simultaneous equations is 
obtained with 'j q. = Q as the independent variable: 



After equations (5) are solved for Li/op, Le/a^ . . . l^s/a^ for 
each value of Q, * 
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2(L , t 

CCq^ cx^p oc^ oc^ j 


( 6 ) 


for the exposed wing. 

The effect of wing-inertia loading acting in opposition to aero- 
dynamic loading has been neglected in equations (s) . The addition of 
a wing-inertia term for a case where the wing is divided into two sections 
is illustrated as follows: 



Axo L-| Lp 1 Lrp / \ 

^ ~ W ^ri®21 + ^2022) 

where 


1 ^ 
KiQ ar 


1 

K2Q ar 




1 


> (T) 


1 


W total model weight 

Wq weight of inboard section of exposed half -wing 

W2 weight of outboard section of exposed half -wing 

Lp total lift of model 
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The inboard movement in the spanwise center of pressure which results 
from the reduction in lift due to aeroelastic wing twist causes a change 
in aerodynamic-center position. The lateral center of pressure of the 
exposed elastic wing can be found from Li/ar, Le/ar, . . , 

solved from the simultaneous equations (5) as follows: 


If the span load is assiimed to move inboard along the half- or 
quarter -chord line, the increment in aerodynamic -center position of the 
complete model due to the inboard movement can be found from the following 
expression: 



^lyi' + L2y2* + + l5y5 ' 

^1 ^ + Li^ + L5 



( 9 ) 


For the tan A use either tan or tan whichever is 


appropriate. 


CONFIDENTIAL 


2k 


CONFIDENTIAL 


NACA RM L52L30 


REFERENCES 


1. Gillis, Clarence L,, Peck, Robert F., and Vitale, A. James: Pre- 

liminary Results From a Free-Flight Investigation at Transonic 
and Supersonic Speeds of the Longitudinal Stability and Control 
Characteristics of an Airplane Configuration With a Thin Straight 
Wing of Aspect Ratio 3 . NACA RM L9K25a, I960. 

2. Glllis, Clarence L., and Vitale, A. James: Wing-On and Wing-Off 

Longitudinal Characteristics of an Airplane Configuration Having 
a Thin Unswept Tapered Wing of Aspect Ratio 3> As Obtained From 
Rocket -Propelled Models at Mach Numbers From 0.8 to l.U. NACA 
RM L50K16, 1951 . 

3- Vitale, A. James, McFall, John C., Jr., and Morrow, John D.: Longi- 

tudinal Stability and Drag Characteristics at Mach Numbers From 0.75 
to 1.5 of an Airplane Configuration Having a 60° Swept Wing of Aspect 
Ratio 2.24 as Obtained From Rocket-Propelled Models. NACA RM L51K06, 

1952. 

4. Chapman, Rowe, Jr., and Morrow, John D. : Longitudinal Stability and 

Drag Characteristics at Mach Numbers From O. 7 O to 1.37 of Rocket- 
Propelled Models Having a Modified Triangular Wing. NACA RM L52A31, 

1952. 

5. McFall, John C., Jr., and Hollinger, James A.: Longitudinal Stability, 

Control Effectiveness, and Drag Characteristics at Transonic Speeds 
of a Rocket-Propelled Model of an Airplane Configuration Having an 
Unswept Tapered Wing of Aspect Ratio 3-0 NACA 65A004.5 Airfoil 
Sections. NACA RM L52L04, 1952. 

6 . Peck, Robert F., and Mitchell, Jesse L. : Rocket-Model Investigation 

of Longitudinal Stability and Drag Characteristics of an Airplane 
Configuration Having a 60° Delta Wing and a High Unswept Horizontal 
Tail. NACA RM L52K04a, 1952. 

7 . Mitchell, Jesse L., and Peck, Robert F. : An NACA Vane-Type Angle-of- 

Attack Indicator for Use at Subsonic and Supersonic Speeds. NACA 
RM L9F28a, 194 9 . 

8. Osborne, Robert S. : A Transonic -Wing Investigation in the Langley 

8 -Foot High-Speed Tunnel at High Subsonic Mach Numbers and at a 
Mach Number of 1.2. Wing-Fuselage Configuration Having a Wing of 
45° Sweepback, Aspect Ratio 4, Taper Ratio 0.6, and NACA 65 AOO 6 
Airfoil Section. NACA RM L50H08, 1950. 


CONFIDENTIAL 


4F 


NACA RM L52L30 


CONFIDENTIAL 


25 


9. Cahill, Jones F., and Gottlieb, Stanley M. : Low-Speed Aerodynamic 

Characteristics of a Series of Swept Wings Having NACA 65A006 
Airfoil Sections (Revised). NACA RM L50F16, 195O. 

10. Parks, James H., and Kehlet, Alan B. : Longitudinal Stability, Trim, 

and Drag Characteristics of a Rocket-Propelled Model of an Air- 
plane Configuration Having a 45° Sweptback Wing and an Unswept 
Horizontal Tail. NACA RM L52F05, 1952. 

11. DeYoung, John, and Harper, Charles W. : Theoretical Symmetric Span 

Loading at Subsonic Speeds for Wings Having Arbitrary Plan Form. 
NACA Rep. 921, 1948. 

12. Loving, Donald L., and Williams, Claude V.: Aerodynamic Loading 

Characteristics of a Wing -Fuselage Combination Having a Wing of 
45® Sweepback Measured in the Langley 8-Foot Transonic Tunnel. 

NACA RM L52B27, 1952. 

13. Diederich, Franklin W., and Foss, Kenneth A.: Charts and Approxi- 

mate Formulas for the Estimation of Aeroelastic Effects on the 
Loading of Swept and Unswept Wings. NACA TN 2608, 1952. 

14. Osborne, Robert S., and Mugler, John P., Jr.: Effects of Wing 

Elasticity on the Aerodynamic Characteristics of a 45° Sweptback- 
Wing— -Fuselage Combination Measured in the Langley 8-Foot Transonic 
Tunnel. NACA RM L52G23, 1952. 

15. Mitchell, Jesse L. : The Static and Dynamic Longitudinal Stability 

Characteristics of Some Supersonic Aircraft Configurations. NACA 
RM L52A10a, 1952. 

Some Recent Data 
NACA RM L50K24, 1951. 


16. Langley Pilotless Aircraft Research Division: 
From Flight Tests of Rocket-Powered Models. 


CONFIDENTIAL 


26 


CONFIDENTIAL 


NACA RM L52L30 


TABLE I.- FUSELAGE NOSE AND TAIL ORDINATES FOR AIRPLANE CONFIGURATION 



in. 

in. 

0 

0.168 

0.060 

.182 

.122 

.210 

.245 

.224 

.480 

.294 

.735 

.350 

1.225 

.462 

2.000 

.639 

2.450 

.735 

4.800 

1.245 

7.350 

1.721 

8.000 

1.849 

9.800 

2.155 

12.250 

2.505 

13 . 125 

2.608 

14.375 

2.747 

14.700 

2.785 

17 . 150 

3.010 

19.600 

3.220 

22.050 

3.385 

24.500 

3.500 
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Figure 1.- General arrangement of model. All dimensions are in inches. 
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(b) Three-quarter front view. 

Figure 2.- Photographs of duralumin- wing model. 
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(a ) Load applied along the 50-percent chord line . 



(h ) Load applied along the 25-percent chord line . 

Figure 3 .- Twist in the free-stream direction per unit load applied at 
various stations along the span of the duralumin wing. 
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(a) 


Before launching 


L-7i^U 


(b) After launching . 

L-77899 


Figure 4.- Model launching. 
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Figure 5.- Static-pressure ratio. 



Figure 6.- Reynolds number of tests, based on mean aerodynamic chords. 
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(a) Steel-Wing model. 


Figiore 7*- Lift data for complete models. 
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(b) Duraliamln wing model. 
Figure 7-- Concluded. 
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M=1.02 M=1.12 M=1.22 M=1.27 

a, deg 

(a) Steel wing. 


Figure 8.- Lift data for exposed wings. 
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(b) Duralumin wing. 
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Figure 8.- Concluded. 
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Figure 9 *- Lift-curve slopes. 




■‘‘^total 



M 


Figure 10.- Rate of change of exposed-wing lift with complete model lift. 
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Figure 11.- Variation of complete model lift-curve slope with wing 

stiffness at M = 1.11. 



-!Slf)rer. 


Figure 12.- Extrapolation "between steel and duralumin exposed-wing lift- 

curve slopes to obtain rigid values. 
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Figure 13 •- Ratio of 
calculated by the 


elastic- to rigid-exposed-wing lift-curve slopes 
method of the appendix for the duralumin wing. 
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Figure 14.- Comparison between experimental and calculated elastic- to 
rigid-lift-c\irve-slope ratios for the duralumin wing. 
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(a) Horizontal-tail deflection, (b) Horizontal-tail deflection, 

-3.0°. 0.14°. 

Figure 15.- Pitching -moment data for the steel-wing model about a center- 
of -gravity position 24.8 percent of the mean aerodynamic chord as 
obtained from the double differentiation method. Flagged symbols 
indicate decreasing Cj^. 
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(a) Horizontal tail deflection, (- 5 ) Horizontal tail deflection, 

-4.60°. -0.72°. 

Figure l6.~ Pitching-moment data for the duralimLn-wing model about a 
center-of-gravlty position 4.60 percent of the mean aerodynamic chord 
as obtained from the double accelerometer method. Flagged symbols 
indicate decreasing Cl. 
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M 

(a) Duralumin-wing model; center-of -gravity position, 4.60 percent mean 

aerodynamic chord. 



.8 .9 1.0 l.l 1.2 1.3 1.4 
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(b) Steel-wing model; center-of -gravity position, 24.8 percent mean 

aerodynamic chord. 

Figure 17.- Period of oscillations. 
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Figure l8.- Variation of aerodynamic -center position with Mach n\miber. 
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(a) Steel-wirig -model time to damp to one-haLf amplitude. 



(b) Duralumin-wing -model time to damp to one-haLf amplitude. 
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(c) Pitching-moment damping derivative 
Figure 19*- Daii5)ing characteristics of the short 


Cmq + Qn^. 

-period oscillations. 


CONFIDENTIAL 


MCA RM 152 L30 


CONFIDENTIAL 


45 




(b) Mach nxjmbers 1.12 and 1.02. 



Figure 20.- Steel- wing -model basic drag data. 
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(b) Mach numbers 1,1k and 1,04. 

Figure 21.- Duralumin-wing -model basic drag data. 
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(c) Mach numbers 0.97 ajcid. 0.94. 



(d) Mach number 0.92. 
Figure 21.- Concluded. 
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Fig\ire 22.- Minimum drag coefficients. 
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(a) Cl ~ 0.50. 
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(b) Cl ^ 0.15. 

Figure 23.- Effect of lift on drag. 


NACA-Langley - 1-28-53 - 325 


COTFIDENTIAL 


SECURITV irsjFORIVIATIOrsi 

CONFIDENTIAL 





Jr* *« 




A 


J 




f 


1 '^ 




» 


1 


« 



I 


* I 

I 




I 


CONFIDENTIAL 




